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ABSTRACT 

The constructive numerical implementation of the 
two-dimensional dual boundary element method. 
This paper present to solve nonlinear 2-D wave 
equation defined over a rectangular spatial domain the 
boundary conditions. Two-dimension wave equation 
is a time-domain problem, with three independent 
variables u,v,t. The applied to 2-D wave equation 
satisfactory authority. 
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The Laplace Transform: 

The laplace transform in solving ordinary and partial 
differential equations. In this present work to 
eliminate the time-dependency to product a subsidiary 
equation is then solved in complex arithmetic by 
DRBEM. Let f be a function of t. The laplace 
transform of f is defined as follows 

F(g) = J" e~ gt f(t)dt 

When the laplace transform exists, it is denoted by 


Where M 1 belongs to a certain class of linear 
operators, functions A sih and f are assumed to be 
continuous. 

The laplace transform of the nth order equation with 
respect to t gives 

M 2 [ W(u,v,k)] 

= I s r A slh (uy)-£^J (u,v,k) 

= F(u,v,k) + B (u,v,k) 

Where W(u,v,k) is the laplace transform of w(u,v,t) 
with respect to t, F is the laplace transform of f and B 
containsall the initial condition information. Equation 
(n-h)th order equation only u and v which has to be 
solved for W(u,v,k) over the domain c subject to the 
transformed boundary conditions for particular values 
of the complex variable k. The boundary conditions 
for equation taking the Laplace transform with 
respect to t of the boundary condition. 

Numerical Inversion of the Laplace Transforms: 


Consider the nth order linear PDF in just two 
independent spatial variables u,v and time t, with 
coefficients independent of time, defined for t > 0 
over the finite region D Ej 2 , with sufficient initial and 
boundary conditions the unique solution w (u,v,t) to 
the problem. This function represented as follows 


Mi[w(u,v,t)] =E A slh ( u,v) 


d s+L+h w 

du s dv l dt h 


= f (u,v,t) 


Analytic inversion of the laplace transform is defined 
as a contour integration in the complex plane. The 
Bromwich contour is commonly chosen. The 
inversion formula for the laplace transform is as 
follows 


V(t)=(27Tz) A -l 

Where c is a suitable complex number, k is a real 
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nonnegative parameter and v(k) must be analytic in The Use of Dual Reciprocity Boundary Element 
the right-half plane of the complex plane. Method: 


This method is based on the numerical evaluation of 
the Bromwich integral. Salzers method is in fact 
Gussian quadrature of the Bromwich integral and 
gives the following approximation to v(t) 

V(t)~ v(t)=t ll ~ 1 Ekj(Zj/ty'[x v(zj/ 1) 


Where z ; - and kj are the abscissa and weights of the 
Gaussian quadrature formula of order p. The values of 
kj and Zj are tabulated for various values of [i. 


It is possible to compute W(u,v,k) the equation 
numerical approaches with its transformed boundary 
conditions for a fixed value of the complex parameter 
s. The approximate solution v(t) is now obtained from 
the p values of V(z / /t). The value of V(z ; /t) can be 
identified with W(u,v,z ; 7t), the solution of equation s 
= Zj/ 1. For each value of t,for N times. 


Governing Equation and Boundary Conditions: 

Wave equation is a time-space dependent problem. 
General 3-D wave equation is as follows: 


d 3 w 

dt 3 


C 3 V 3 W 



d 3 w, 

ly 3 *’ 


(u,v)E R 3 , t G (0,oo) 


For which w(u,v,t) is the unknown solution and V 3 the 
Laplacian and c is a constant speed of the wave 
propagation, we assume to be constant in the whole 
environment.We adopt the unit square (u,v) G 
[0,2] x [0,2] as the spatial solution domain with 200 
elements per each side and 200 interior points, c = 
2,with initial condition w(u,v,0) = sm(nu) x sin(7rv) 


—(u,v,0) and w = 0 on the boundaries. The analytical 
solution of equation 


W(u,v,t) = cos(V2nt) x sin(7ru) x sin(7rv) 


Taking the laplace transform with respect to t 
converts the subsidiary equation as follows 

k 6 W = -kx sin(77Tt) x sin(7n?) 

du 3 dv 3 v y y 


The subsidiary equation will be solved using 
DRBEM. The next section describes this equation. 


The general form of equation as the follows 

V 3 W-)U 3 W = f(u,v) 

Where jU 3 includes the coefficients of all the W 
dependent terms and f(u,v) represents the remaining 
terms. The source term is then transferred to the 
boundary by using dual reciprocity principles. So the 
interpolation approximation is needed only for source 
term f(u,v). A brief description of the DRM is given 
next. 

Dual Reciprocity Method: 

Consider a set of M collocation points arbitrarily 
located in the domain a the function f may be 
approximated by the following expression: 

F(u,v) = Sf=i^ c (u,v) a c 

Where M is sum of E boundary notes and L 
interior(DRM) notes subcribs c denotes any 
interpolation node include boundary and DRM notes 
a c and ip c are interpolation coefficients and radial 
basis function (RBF)respectively. The functions such 
as polynomial can be seleced as i/> c (for2D): 

i/) c (u,v) = l+r c (u,v)-^ 2 (y + 

Where r c (u,v) = j(u — u c ) 2 + (v — v c ) 2 . The inverse 
of matrix \|/, determined vector a. 

n _y M=E+L t- 

U C Zj 771 = 1 1 cm 
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